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Abstract 

It is shown that fractional derivatives of the (integrated) invariant measure of the Feigen- 
baum map at the onset of chaos have power-law tails in their cumulative distributions, whose 
exponents can be related to the spectrum of singularities f(a). This is a new way of char- 
acterizing multifractality in dynamical systems, so far applied only to multifractal random 
functions (Frisch and Matsumoto (J. Stat. Phys. 108:1181, 2002)). The relation between 
the thermodynamic approach (Vul, Sinai and Khanin (Russian Math. Surveys 39:1, 1984)) 
and that based on singularities of the invariant measures is also examined. The theory for 
fractional derivatives is developed from a heuristic point view and tested by very accurate 
simulations. 

Keywords: chaotic dynamics, multifractals, thermodynamic formalism. 

1 Introduction 

Recently a new method for analyzing multifractal functions u(x) was introduced pQ. It exploits 
the fact that the fractional derivative of order a (denoted here by D a ) of u(x) has, for a suitable 
range of a, a power-law tail in its cumulative probability 

Prob{|L> a u| > £} oc«T p *, f -» oo. (1) 

The exponent p*(a) is the unique solution of the equation 

C P * = P*(a)a, (2) 

where £ p is the scaling exponent associated to the behavior at small separations I of the structure 
function of order p, i.e. (\u(x + l) — u(x)\ p ) ~ It was also shown that the actual observability 
of the power-law tail when multifractality is restricted to a finite range of scales is controlled by 
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how much ( p departs from linear dependence on p. The larger this departure the easier it is to 
observe multifractality. 

So far the theory of such power-law tails has been developed only for synthetic random 
functions, in particular the random multiplicative process [2] for which Kesten-type maps [3] 
and large deviations theory can be used. 

It is our purpose here to test the fractional derivative method for invariant measures of 
dissipative dynamical systems, in particular for the Feigenbaum invariant measure which appears 
at the accumulation point of the period doubling cascade where the orbit has period 2°° |IJ[3]. 
Its multifractality was proven rigorously in Ref. |Bj using a thermodynamic formalism. For the 
Feigenbaum measure all scaling exponents can be determined with arbitrary accuracy. 

There is an important difference in the way one processes functions and invariant measures 
to determine their multifractal properties and in particular the spectrum of singularities, usually 
denoted D(h) for functions [7] and f(ct) for measures [HIE]. For a function u(x) one uses the 
moments or the PDFs of the increments u(x + I) — u(x) to determine the scaling exponents, 
whereas for an invariant measure fio one works with integrals over intervals or boxes of different 
sizes. In the one-dimensional case the two approaches become equivalent by introducing the 
cumulative distribution function 



Hence we shall apply the fractional derivative method to the integral of the invariant measure. 

The organization of the paper is the following. Section |2 is devoted to the thermodynamic 
formalism for the Feigenbaum attr actor. In Section l2.1( we recall the method used in Ref. [Sj. 
In Section 12.21 we show how this formalism, based on the study of the geometrical properties 
of the attractor, is actually connected to the standard multifractal formalism which focusses on 
the statistical properties of the invariant measure [HJ|U]. To the best of our knowledge the exact 
relation between the two formalisms is discussed here for the first time. Then, in Section [2.31 we 
calculate numerically the free energy and accordingly the scaling exponents ( p for the integral of 
the invariant measure; this is done by a very accurate transfer-matrix-based method. Fractional 
derivatives are discussed in Section 01 In Section 13.11 we briefly recall the phenomenology of 
power-law tails in the distribution of fractional derivatives and the limits on observability. The 
fractional derivative analysis of the Feigenbaum measure is presented in Section 13.21 Concluding 
remarks are made in Sectional 

2 Thermodynamic Formalism for the Feigenbaum Attractor 
2.1 Thermodynamic formalism 

In this section we give a brief description of the thermodynamic formalism for the invariant 
measure of the Feigenbaum map (see Ref. jHJ for the mathematical details) and show how one 
can use it in order to study the multifractal properties of the Holder exponents. 

By Feigenbaum attractor we understand the attractor of the one-dimensional mapping g : 
[0, 1] — > [0, 1] , where g{x) is the solution of the Feigenbaum-Cvitanovic doubling equation: 




(3) 




(4) 
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Equation Q is known to have the unique solution in the class of smooth unimodal maps (that 
is, maps having one critical point) with a non-degenerate maximum. This solution is called the 
Feigenbaum map. It is holomorphic in some complex neighborhood of [—1, 1] and the first few 
terms in the power series expansion are |lf)| 

g{x) = 1 - 1.5276 . . . x 2 + 0.1048 . . . x 4 + 0.0267 . . . x 6 - 0.003527 . . . x 8 + . . . . (5) 

The value of the universal constant 7 which is the inverse of the Feigenbaum scaling constant a is 
approximately equal to 0.3995. An attractor A for the map g can be constructed in the following 
way. For each n > 1 define a collection of intervals of n— th level: 

A< n W°(A£ n) )=0O0o...o 5 (A< n) ) (l<i<2«-l). (6) 

i 

The following properties of the intervals A^ are easy consequences of the doubling equation 
®: (a) Intervals Af°,0 < i < 2 n - 1 are pairwise disjoint, (b) g^A^ C A[, n) . (c) Each 
interval of n— th level A^ contains exactly two intervals of (n + 1)— th level, A 4 - n+1 ^ and A^?J„ . 

(d) 27 2n < lA^I < 27", where | • | denotes the length of the interval. The first three levels of the 
intervals are shown in Fig. ^ 

The properties above imply that it is natural to use a dyadic representation for the intervals 
A 4 - n \ Let i = Y^Jj=l € j 2*, where €j = 0, 1. Then we can use a sequence (eo, . . . , e n ~i) as a symbolic 
coding for intervals A^: A^ = A^„ en _ 1 . Now we can define the Feigenbaum attractor 

^nuV- ( 7 ) 

n>l i=0 



3 



The set A is isomorphic to the set of all infinite dyadic sequences (eo, e±, . . . , e n , . . .). Such se- 
quences can be considered as a symbolic coordinate system on A. In this new coordinate system 
the map g acts as the dyadic addition of the sequence (1,0,... ,0,...). Notice that topologically 
A is a Cantor set. It is easy to see that A is indeed an attractor for all but countably many initial 
points a: €[-1,1]: dist(c/( m ) (x), A) -> as m — > oo. The exceptional set of initial points consists 

of all unstable periodic orbits and their preimages. 

(n) 

As we have seen above, all intervals AJ have exponentially small lengths but the exponent 

varies from 2 7 n to 7 2n . Notice that exponents ln|A- n ^| give all possible scalings of the fractal 

set A. The basic ingredient which is needed for the multifractal analysis is the control over the 

(n) 

spectrum of possible scalings corresponding to exponents, i.e. hi|A^ |. Such control can be 
achieved with the help of the thermodynamic formalism. The thermodynamic formalism which 
was constructed in Ref. [0] is based on the Gibbsian description for the lengths of the intervals 
A. n \ It is shown in Ref. that there exists a function U{e^\ , . . . , e^ n \ . . .) (thermodynamic 
potential) which is defined on all infinite dyadic sequences such that: 
I. There exists a constant C > for which 



1 

C 



(1) f (2) 

C 1 



(n) 



II. For any two dyadic sequences 



(n) -(n+1) -(n+2) 



■ )<C. 



(8) 



In) ~(n+l) ~(n+2) 
, c , c , c , 



which coincide on the first n positions 
|£7(eW , , . . . , e («) , e -("+D , e ("+ 2 ) ,...)- U(e^ , e^ 2 ) , . . . , e W , e , e ( n+2 ) ,...)!< C(2 7 ) n . (9) 
III. For any Ago J ... )£n _ 1 with eo = 1 



exp[-C7(2 7 )3] < 



IA 



(n) 



£Oi-i e n-l I 



IA 



(n-1) i 
eo,---,e n -2 1 



exp(-C/(e n _ 1 , . . . , ei, 1, 0, . . . , 0, . . .)) < exp[C(2 7 ) a ] . (10) 



It immediately follows from (jlUj) that for Ci = exp 



IA 



c 



l-(2 7 )3 



< 



(n) | 
l,ei,...,en-i I 



exp 



YT s =l U{ts, e a -i, . . • , ei, 1, 0, . . . , 0, . . .) 



The condition eo = 1 corresponding to intervals A^ with odd i's plays only a technical role and 
it is not essential for our further analysis since the odd intervals contain information about the 
lengths of the even ones. Indeed, it is very easy to see that for every odd i the intervals A^ and 

have lengths of the same order. 
We next introduce a parameter (5 (inverse temperature) and define the partition function 



<Ci. 
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Z n(P) = X] eXp 
e n -l,.--,<;o=0,l 



n-1 



/3^C/(e s ,e s _i,...,ei,l,0,...,0,...) 



s=0 



(12) 
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and the free energy 



FM = hMPl, m= lim FM . 



n 



It immediately follows from (|12|) and (|13|) that 



2 n -l 



Wi/3 



exp[nF (/?)]. 



(13) 



(14) 



i=0 



In the thermodynamic limit n — > oo the probability distributions 

n-l 



1 



exp 



/3^J7(e s ,e s _i,...,ei,l,0, 



,0,. 



s=0 



(15) 



tend to a limiting distribution ug which can be considered as a Gibbs measure with the potential 
U, inverse temperature j3 and the boundary condition eo = 1, e_j = 0, i £ N. This Gibbs 
distribution generates the probability measure on A + = A n A^ which is the part of the whole 
attractor A corresponding to intervals A^ with odd numbers i. We shall denote this Gibbs 
measure on A + by ug. Notice that /3 = corresponds to a unique invariant measure and (3 = 1 
gives a conditional distribution corresponding to Lebesgue measure on [—1,1]. The free energy 
F(f3) contains all information about the multifractal properties of the Feigenbaum attractor. 
Notice that the thermodynamic formalism leads to one-dimensional statistical mechanics with 
exponential decay of interactions and hence without phase transitions. This implies that F{(3) 
is a smooth function. In fact it is holomorphic in some complex neighborhood of the real axis. 
Denote 



n-l 



H(e n - 1 ,...,ei) = ^C/(e s ,e s _i,..., ei,l,0,...) 



(16) 



Using relations 



s=0 



77m Y] -#(e n _i,...,ei)exp[/3#(e 

e n _l,...,e =0,l 



n— 1 1 



In I A 



(n) 

l,ei,...,e n _l I 



tl 



(n) 



+ 0(i 

n 



(17) 



and 



nKy > e„_i,...,e =0,l v 



e„_i,...,eo=0,l 



ei) exp[/?i?(e n _i,...,ei)] 



In I A 



(n) 

l.ei ,...,e r , 
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,ei)exp \j3H(en-i, 
» i \ 



In I A 



l,ei,...,e„_i I 



(n) 



+ o(- 

.n. 



(18) 
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we conclude that F{(3) is a monotone decreasing convex function. We shall also use the spectral 
representation for the free energy. Consider the transfer-matrix operator C{(3): 

C(p)h(eW , , . . . , ,...)= £ exp\J3U(e^ , e« , e< 2 > , . . . , , . . .)}h(e^ , e« , , . . . , , . . .) 

e(0)=0,l 

(19) 

Since £((3) is a positive linear operator, its largest eigenvalue A(/3) is strictly positive and simple. 
It is easy to see that 

F(f3) = \n\(f3). (20) 

For an arbitrary point x E A+, denote by A^ n \x) the interval of the nth level which contains 
x. It follows from (|17|) that for points x which are typical with respect to U/j (that is corresponding 
to a set of full /^-measure.) 

\£S n) [x)\ ~ exp[F'(/3)n] . (21) 
More precisely, for /jg- almost all x G A 

In IA(")(V)I 

Jim ln|A m = F'(P) . (22) 

n^oo 77, 

We next find the total number N n {(5) of the intervals of n— th level whose length is of the order 
|A( n )(/3)| = exp[F'(/?)n]. We have 

N n (P)\A^((3)f ~ JVnOSJexp^GSjn] ~ Z n (/3) ~ exp[F(/J)n] (23) 

which gives 

iV n (/?) ~ exp[(F(/3) - PF'((3))n] . (24) 

Using (|24[) we can find the Hausdorff dimension dn(P) of the set of points x 6 A + which are 
typical with respect to the measure Since 

N n {p)\£sF>{p)\ dB W ~ exp[(F(/3) - (3F'(P) + efo(/3)F'(/?))n] (25) 

we conclude that 

- /?F'(/3) + d H ((3)F'{(3) = (26) 

which immediately implies 

*hW- ^ -^-J^y (27) 

The Hausdorff dimension dn(A) of the whole attractor A is equal to the maximum of dn(f3) 
over all (3 € K. Let /3o be the unique solution of the equation = 0. It is easy to see that 

d H (A) = d H ((3 ) =Po- 

We next discuss multifractal properties associated with the Holder exponents. Consider the 
integral u(x) of the invariant measure uo, defined by (j3J), which is plotted in Fig. [21 The attractor 
being topologically a Cantor set, u{x) is a variant of the Devil's staircase (see Ref. Section 
8.2). 
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-1 -0.5 0.5 1 



x 

Figure 2: The integral u(x) of the Feigenbaum invariant measure calculated with 2 15 bins of 
uniform length A = 2 -14 in [—1, 1]. Inset: The invariant measure smoothed over the distance A 
calculated as a frequency histogram. 



To find its spectrum of Holder exponents, notice that for every interval A^ the increase of 
u(x) along the interval is equal to 2~ n . Hence A. n ^ corresponds to a Holder exponent h n (i) = a 

N I 2" 



where |A- n ^| Q ~ J^. This implies 



|Af°| ~ 2-S. (28) 

Using Q21|). (|2*H)) we conclude that the Holder exponent a corresponds to an inverse temperature 
/3(a) such that 

F\f3{a)) = -— . (29) 
a 

This gives 

0(a) = (FT 1 ("^) , (30) 

where (F')~ (y) is the inverse function to F'(x). We can now find the Hausdorff dimension djj(a) 
of the set of points x £ A for which the Holder exponent of u(x) is equal to a: 

d H {a) = d H ((5(a)) = (5(a) - . (31) 

Notice that the analysis presented above can be made completely rigorous (see, for example, 
Refs. 0H2]). 



2.2 Connection between thermodynamic formalism and standard multifractal 
analysis 

It is quite interesting to compare the multifractal analysis which we presented above with the one 
introduced in Ref. 0. Although we restrict ourselves here to the case of the Feigenbaum attractor, 
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the result presented below holds in a much more general setting. Basically our analysis is valid 
whenever the system under consideration can be described with the help of the thermodynamic 
formalism. 

The basic object for our analysis is the Feigenbaum attractor itself and the method is based 
on the construction of the thermodynamic formalism for the lengths of the elements of dynamic 

(n) 

partitions A^ . The thermodynamic formalism uses considerable dynamical information about 
the map g. In contrast, the analysis in Ref. [S] is carried out for fractal measures and does not 
directly use the dynamical information about the system. In the period-doubling setting the 
fractal measure is fiQ. It is the unique invariant measure for g acting on A (see Ref. j^]). It also 
can be considered as a physical or Sinai-Ruelle-Bowen (SRB) measure on [—1,1]. This means 
that under dynamics given by the map g any initial absolutely continuous distribution //on [—1,1] 
converges to /iq: g^ n \ji) — » hq as n — > oo. The multifractal analysis in Ref. f9j is based on a 
function D p which can be defined in the following way. Consider a partition of the interval [—1,1] 
into subintervals Ij, 1 < j < N of length I. Then 

' " 1 F((3(a)) . (32) 




Inl I ln2 

It follows from (j3H) that 

TV 

J>o(W ~ . (33) 

i=i 

Another characteristic of a multifractal measure is given by its spectrum of dimensions f(a) 
which is just the Legendre transform of D p (p — 1): 

/(a) = inf [op - D p (p - 1)] . (34) 
p 

The dual Legendre relation allows one to find D p (p — 1) from f(a): 

D p (p - 1) = inf [ap - f(a)] . (35) 

a 

We next find a correspondence between the pair {D p (p— 1), f(a)) and the pair of thermodynamic 
functions (F(0), cZ#(a)). We shall show that 

D p (p-1) = -F' 1 (pln2), f(a) = d H (a), (36) 

where i 7-1 is an inverse function to the free energy F(f3). To derive the first relation we consider 

(n) (n) 

the dynamical partition , < i < 2 n — 1 and assume that n ^> 1 but / <C minj \A\ \. For 
each A^ define 

M(A^)= MW- (37) 

Notice that the asymptotic behavior of M(A^) depends only on asymptotic scalings of smaller 
elements of the dynamical partitions Ai n+m ' ) inside A^ n \ The thermodynamic formalism con- 
structed above implies that asymptotically those scalings are completely determined by the po- 
tential U and hence they do not depend on i. Rescaling the invariant measure inside Af 1 ^ by a 
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factor 2 n we conclude that 



„'.r.^A< n ) 



-JV 



D g (l-P) 



(38) 



where JVj is a total number of the intervals L inside A- n \ Taking the sum over i and using (JT 
we have 



E(w>(W = E M ( A H ~ 2^ ^ "~2»/< 

j=l 0<«<2"-l 0<t<2"-l 0<i<2™-l 



J_ y- ^,(i-p) „ J_ 



I A 



(n),\ d p(!-p) 



±. N D p{ l- P ) eMnF{Dp{1 _ p))l 



This together with (|33j) immediately gives 

exp[F(D p (l-p))]=2P 



(39) 



(40) 



which implies the first relation in (|3fij). We next show that the second relation holds. Using (j34j) 
we have 



/(a) = inf [ap - D p (p - 1)] = inf [ap - (-F -1 (pln2)] = inf 

"In 2 



inf 



a 



F(J3) + (3 



p 
a 



a i 



inf 

In 2 /3 



a 



-P + F(f3) 



(41) 



/3 Lin 2 

It is easy to see that the extremum in (jiTjl corresponds to /3(a) = {F')" 1 ( — ^r) which implies 

(42) 



o 



/(a) = /3(a) + _F(/3( a )) = d H (a) . 

Finally we express the scaling exponents £ p for the structure functions through the thermo- 
dynamic characteristics. The exponent Q p is defined by the scaling relation 

(43) 



in terms of the integral u{x) of the invariant measure. Let Ij, l<j<N=jbe& partition of 



(\u(x + I) - u(x)\ p ) ~ 
! in 

[—1,1] into intervals of length I. Then 

1 N 

(| u(x + /) _ u(a .)|P) = ^(^(/.Or - ZiV^C 1 ^) ~ i*+A.(P-i) (44) 

which together with (|3f))) gives 

C p = l + Z? p (p-l) = l-F- 1 (pln2). (45) 
Using (|35|) one can also write £ p in the following form: 

C p = 1 + inf [ap — f(a)] = 1 + inf [ap — (i#(a)] . (46) 
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At a = a* the inflmum is attained, we can then write the relation between p and a* as 

F(/3K)) 

P = atf(")U=a* = — • ^ 

The scaling exponents for the structure functions are hence obtained as 

C P = l- /?(«*) • (48) 



We now turn to concrete calculations for the scheme presented above. It is easy to see that 
all the thermodynamic functions can be effectively approximated numerically. The first step is 
to find approximations for the thermodynamic potential U(e^\ e^ 2 \ . . . , e- n \ . . .). We shall use 
Markov approximations U k (e^ , e^ 2 \ . . . , e^) which are defined by the following formula: 

iaW I 
U k ^\e^\...,e^) = lim In ' y^'" (49) 

1 l,0,...,0,e( fe ),e( fc - 1 ),...,e( 2 ) 1 



It was shown in jHj that the limit in (|49|1 exists and 

\U(e^\e^\...,e( k \e^ 1 \...)-U k (e^,e^,...,e^)\<C(2 1 )K (50) 

Using the approximate potential U k we can construct a Markov approximation for the transfer- 
matrix operator C{(3). Namely, we define a finite dimensional linear operator C k {(3): 

£ fc (#M* ( V a) ,...,e W )= E exp^x^^^ 

e(°)=0,l 

(51) 

In the matrix representation the operator C k {f3) corresponds to a certain 2 k x 2 k matrix. Denote 
by \k(P) its largest eigenvalue. Then 

F k ((3) = In \ k (p) (52) 

is a natural approximation for the free energy F((3). It follows from (|50|) that F k {(5) converges 
to F(f3) exponentially fast in C°° topology. Using F k {j3) we can effectively approximate all 
the multifractal functions which we discussed above. The corresponding numerical results are 
presented in the next section. 

2.3 Numerical calculation of the free energy and the scaling exponents 

Here we show how to construct the transfer-matrix operator C k starting from k = to general 
k. For k = 0, the matrix operator Co is just a scalar. Denoting A^... i£n by A^ n ^(ei, . . . ,e n ) for 
clarity, let us consider 

n-2 

exp,^'"')] = lim ^M^M, (53) 
| A^- 1 ) (1,0,..., 0)| 

n-2 
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whose analytical expression is easy to calculate. The 0-th order approximation of the free energy 
is given by 



Fo(^)=ln[ e ^(°) + e ^W]. 
For e(°) = 0, since A { ™ ] = [-a~ n , a~ n ] and = [g(a~ n ), 1], we have 
|A^>(1,0, ... ,0)| = l-g(a- n ) = ca~ 2n + 0{a~ 4n ), 



(54) 



(55) 



where we use the expansion © of g(x). Thus c = 1.527.... Therefore, the corresponding 
component of Cq(0 = 1) is 

\ca~ 2n + 0(a" 4n )| 



=A(o) = lim 



a 



(56) 



n^oo | ca -2(n-l) + 0( Q -4(n-l)) | 

Before considering the second term in the argument of the logarithm in (|54[) , we recall the relation 

g^){a- n x) = {-l) n a- n g{x), (57) 

which played a primary role in the proof of Theorem 4.1. of Ref. 0. Since AW(1,0,...,0,1) = 
A^ 2 „_ l5 in view of (|57|) . we have 

|AW(1,0,... ,0,1)1 = \g^\ti)-gV+^\a")\ 

= \ca-^ n -^{\g{a- x )] 2 - 1} + 0(a~< n -V)\. 

Thus for e (°) = 1, © can be rewritten as 



(58) 



1 - [g(a 



-1m2 



|ca-2(n-l) +0( Q -4(n-l))| 

We now arrive at the expression for the 0-th order approximation of the free energy: 

F (P) = \n[a- 2 P + {l-[g(a- 1 )] 2 n 



(59) 



(60) 



Next, consider the k = 1 approximation to the transfer matrix C\{j3). From (|51|). it can be 
written in standard matrix notation as 

e /3(7 2 (0,0) e /3U 2 (l,0) 
It then follows that the free energy is given by 



(61) 



In 



e 

In 2. 



/3(7 2 (0,0) + e /3(7 2 (l,l) + ^/( e /3l/ 2 (0,0) _ e /3l/ 2 (l,l))2 + 4 e /3[/ 2 (l,0) e /3[/ 2 (0,l) 



(62) 



The four components of the transfer matrix require the evaluation of suitable exponential terms, 
expressible by (|49|). from 



exp[£/ 2 (e (0) ,e (1) )] = lim 



n— 3 

|AW(l,0^0,eW,e(°))| 
|A(^ 1 )(l,p,...,q ) e( 1 )] 

n— 3 



(63) 
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Each of these terms is calculated in the same manner as for the k = case: 

^2 (0,0) 



v |A(")(1,0,...,0)| 2 
hm — -. — -T- — = a , 



^2(1,0) 



n^oo | A^- 1 ) (1,0,..., 0)| 

|AW(1,0,... ,0,1)| 



,tf2(l,l) - 



lim 

n^oo 0, . . . , 0)| 

1 " [5(«- 1 )] 2 , 

|A(")(1,0,..., 0,1,0)| 

n-oo ^("-^(l.O, ...,0,1)| 

|b(0)] 2 -[g(q- 2 )] 2 | Hg(gif 
|b(0)]2-[ 5 (a-i)]2| i_[ 5 ( a -i)]2' 

|A(")(1,0,..., 0,1,1)| 

n^oo jA^-^^O, ...,0,1)| 
| b (3) (0)] 2_ b (3) (a - 2)] 2| 



l\i2 



\g(0)Y ~ MO] 



(64) 



(65) 



(66) 



\[g(0W - [g(a^W\ ' 
For = 2, we just write down the transfer matrix £2(0) '■ 

I e f}U 3 (p,0,0) e /3C/ 3 (l,0,0) q 

e /3C/ 3 (0,0,l) e /3C/ 3 (l,0,l) 

e /3C/ 3 (0,l,0) e /3E/ 3 (l,l,0) 

o e /3c/ 3 (o,i,i) e /3f/ 3 (i,i,i) y 



(67) 



(68) 



Now we are in a position to calculate the transfer matrix £k((3) for general k. Let us consider 
the component for (3 = 1: 



exp[£4 +1 ( e W e 



(0) ,(i) 



lim 



n-k-2 

|A(™)(l,oT^To,e (fc) ,---,e (1) ,e (0) )| 



n-voo |A(^- 1 )(1,0, . . . ,0,eW, . . . ^C 1 ))! 

n-k-2 
n-k-2 

|g(AW(C^,eW,...,e( 1 ),e(°)))| 
™ |#(A( n_1 )(0, . . . ,0,eW, . . . ,e^))\ 

n-k-2 



(69) 



Hence it is enough to calculate |A^(0, . . . , 0, e^ k \ . . . , e^, e^)|. For this we use the following 
relation: 

|A<»>(0^ £ W ..., e <V 0) )l 

n-k-2 

_ |^( e ( fc )2 n - fe - 1 + ...+e( 1 )2 n - 2 +e(°)2 n - 1 )^g-j _ g(e( k '>2 n - k - 1 +...+e( 1 '>2 n - 2 +e < - '>2 n - 1 ) (q,-™) | 

= a-("- fc )|A( fc+1 )(e( fc ),...,e( 1 ),e(°))|. ' (70) 
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By (|70[) the numerator in (J69|) is given by 



IA 



(n) 



|{1 - ca- 2 ^- 1 )^)^)] 2 } - {1 - ca- 2 ^ k -^[g^(a- {k+1 ^} 2 } + 0{a~ A ^)\} 



= ca- 2i ~ n ~ k ^\[g^m 2 - b (j) (a~ (fc+1) )] 2 | + O^" 4 ^). 
Therefore the component (|6l?]) is expressed as 

V>(o)] 2 -b&V- (fc+1) )] 2 



(71) 



exp[C/ fc+1 (eW,e 



(0) ,(i) 



[ 5 0-')(0)]2-[ 5 0-')( a -fe)]2 



(72) 



where j 



e (k) + e (fc-l) 2 + . . . + e ( 1 )2 fe - 1 + e(°)2 fc and j' = + e^ 1 ^ + . . . + e^2 k ~ 1 . In other 

„(*+*) 



, ■ - • (k-\-l) 

words, if the coordinates of the end points of the intervals are given by 



and A^ 



3 



b - t ], the component (|69j) can be rewritten as 



exp[C4+i(e 



(0) ,(1), 



ffc+i) 



(fc+i) 



,(fc) 



(73) 



If the transfer matrix is written as an ordinary 2 k x 2 fc matrix, it is easily found that the row and 
columns indices are given by (row, column) = (p,q), where 



(row) p 
(column) q 



l + e (*) +e (*-i) 2 + ... + £ (i) 2 



k-l 



1 + e (fc-i) + e (fc-2) 2 + . . . + e 



(O)oJfe-l 



(74) 
(75) 



For numerical calculation of the functions Fk(f3) and £ p associated to u(x), we use the ex- 
pansion (JSJ) for <7(x), as given in Ref. jlUj . We include terms up to x 32 . Numerical calculation is 
done with standard double precision (15 significant digits). For obtaining the largest eigenvalue, 
we use the power method for matrices |13j . The multiplication of the matrix is stopped when 
the relative error of the most dominant eigenvalue becomes less than 1CP 13 , thus giving 1CP 13 
absolute error on the F{(3) function. An alternative approach, also based on the thermodynamic 
formalism and yielding 10-digit accuracy, may be found in Ref. [12]. Approximate free energies 
Fk(/3), with k up to 16, calculated by the transfer-matrix method are shown in Fig. [3] We note 
that the k = approximation (|60|) already gives a reasonable estimate. The discrepancy of the 
free energy between various orders of approximations is visible at large (5. However the (3 > 
region is irrelevant as far as ( p for positive p is concerned (see Eq. (|48|0 . The corresponding 
Cp (p — 0) are calculated from the Fk(f3) by for different values of k; the results, which hardly 
depend on k, are shown in Fig. [If a). We also determined the structure functions of u(x) with 
2 17 uniform bins in [—1,1]; they are plotted in Fig. EJb). The exponents Cp are then obtained 
by a least square fit of the structure functions over the range 2~ 16 ~ 1.5 x 10~ 6 < r < 1. With 
this number of bins, the quality of the fit begins to somewhat deteriorate beyond p = 4, but 
otherwise there is rather good agreement between the two methods of determining £ p . Note that 
the a y- intercept" of the graph of £ p , namely (o, which is the codimension of the support of the 
invariant measure no(dx), is positive and its numerical value is slightly under one half 6i. This 
will be important in the sequel. 
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20 




J8 



Figure 3: The k-th order approximation to the free energy Fk{(3). Inset: Enlargement of the 
region 6 < f3 < 10; when increasing k, alternate convergence is observed. 



2 - 





(a) 


















































k=i 
k= 

k=: 


) 

L - 
I 
















k=4 

A-8 - 
k=12 - 
k=16 
asured ° 














me 



5 
P 



10 




Figure 4: (a) Scaling exponents £p of structure functions obtained by two methods. Open circles: 
data obtained by a least-square fit of the slopes of the directly measured structure functions 
shown in (b). Lines: data calculated from the free energy Fk{(3) using (|45|) (curves for different 
k are essentially indistinguishable). 
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3 Fractional derivatives for the Feigenbaum attractor 



3.1 Phenomenology for multifractality and fractional derivatives 

In this section we briefly recall the phenomenological approach to multifractality via fractional 
derivatives and adapt it to a multifractal measure. We therefore work, not with the measure 
Ho(dx) itself, but with its integral u[x). Singularity exponents a may be viewed as local Holder 
exponents of u(x), i.e., \u(x + I) — u(x)\ oc \l\ a for (I — ► 0). We turn to fractional derivatives 
of order a defined, as in Ref. pQ, as the multiplication in the Fourier space by u(k) by \k\ a (see 
Ref. for precise definition). An isolated non-oscillatory singularity with exponent a at a point 
x implies 

\D a u(y)\ ~ \y - x\ a ~ a (y -> x). (76) 

If a — a < 0, as we shall assume hereafter, the exponent is negative, the fractional derivative 
can become arbitrarily large and thus contributes to the tail-behavior of the probability. A key 
assumption in the phenomenology is that this argument can be carried over to non-isolated multi- 
fractal singularities, provided we take all types of singularities into account. For the Feigenbaum 
invariant measure, we know the Hausdorff dimension djj(a) of the set of points having a singu- 
larity with exponent a. Assuming that we can also use dnioi) as a covering dimension, we can 
express the probability to have a singularity of exponent a contributing a fractional derivative of 
order a which exceeds (in absolute value) a given large value £, that is we require 

i 

\y — x\ < £ a ~ a . (77) 
In terms of the codimension of the set I a , the probability to satisfy (f77j) is written as 

Prob{|L> a ii| > oc \y - x \ d - dH{a) cx f~^5T\ (78) 

Here d is the spatial dimension (d = 1). Taking now into account the singularities with all possible 
exponents a, the tail of the cumulative probability of the fractional derivative of order a is given, 
to the leading order, by the following power law 

Prob{|£> a u| > £} oc C p * , £ ^ co, (79) 

. d-d H (a) 

p* = mf . (80) 

a<a a — a 

An easy calculation shows that a* corresponding to the infimum in l|8fljl satisfies 

d H (a*) - d 

a+ = a*(a) = a-\ -— — r — , (81) 

Off (a*) 

which immediately gives = d' H (a±). On the other hand, we know that 

( p = inf (pa + d- d h (a)). (82) 

a 

Here the infimum is given by an a satisfying the very same relation = d' H (a ic ). Hence, 

C P * = + d- d H (a±). (83) 
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Using ([STjl. we get 



C P * = P* 



a + 



d H {a*) - d 



) 



+ d- d H {a*) = P*a, 



(84) 



where the second relation follows from p± = d' H (a*). The geometrical interpretation of this 
equation is that the (negative) exponent of the power-law tail for the fractional derivative of 
order a is the p-value of the intersection of the graph of C p and of a straight line of slope a 
through the origin. 

As shown in Ref. pQ, in the presence of the finite range of scaling, the power-law tail (|79|) 
emerges only if the multifractality is sufficiently strong. This strength is given by the multifrac- 
tality parameter C{a) , a measure of how strongly the data depart from being self similar (which 
would imply Q p oc p) : 



where a* = d(^ p /dp\ p=Pit . It was shown that observability of the power-law requires a sufficiently 
large value for the product nC(a), where n is the number of octaves over which the data present 
multifractal scaling. In practice it was found in Ref. that 



For example, fully-developed turbulence velocity data JH] have typical C values of the order 
of 1/30, thereby requiring a monstrous inertial range of about 300 octaves for observability of 
power-law tails. As we shall see, the situation is much more favorable for the Feigenbaum invariant 
measure. 

Before turning to numerical questions, we comment on an issue raised by an anonymous referee 
who worried about the nonlocal character of the fractional derivative and wrote in essence that our 
approach makes sense, strictly speaking, only for (statistically) translationally invariant in space 
systems: otherwise, if the system consists of components whose "fractal properties" are rather 
different the results will be smeared out. Our feeling about such matters is summarized as follows. 
First one can observe that, of course, the attractor for the Feigenbaum map is not homogeneous 
(translation-invariant) but after zooming in it becomes increasingly so; the fractional derivative 
is not a local operator but the tail of its PDF is likely to be dominated by strongly localized 
events. Second, a more technical observation. The idea of the multi-fractal analysis is based on 
the fact that the dynamics of a system determines a variety of scales. It is important that these 
scales do not depend on a particular place in the phase space. On the contrary, they are present 
and "interact" with each other everywhere. In the case of the Feigenbaum attractor the scales 
depend on a symbolic location in a system of partitions. In physical systems, like homogeneous 
turbulence, such partitions are difficult to define rigorously. However, the invariance with respect 
to the space coordinate is still present and forms a basis for applicability of the multifractal 
calculus. 

3.2 Numerical analysis of fractional derivatives 

The phenomenological arguments presented in the previous section suggest that we should find 
power-law tails in the cumulative probability for fractional derivatives of u{x) for suitable orders a. 




(85) 



nC{a) > 10. 
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Table 1: For the Feigenbaum invariant measure we show the scaling exponents £ p , the corre- 
sponding inverse temperature /?'s, the multifractality parameter C and the number of scaling 
octaves needed. 



p 


C P 


P 


C 


number of octaves needed 


1.0 


1.0 


0.0 


0.48 


21 


1.5 


1.2540292658 


-2.54029265895 


0.34 


29 


2.0 


1.4985620106 


-4.98562010659 


0.28 


36 


2.5 


1.7344372955 


-7.34437295506 


0.24 


42 


3.0 


1.9625763533 


-9.62576353333 


0.22 


45 



Inspection of Fig. 0] indicates that a should be between the minimum slope of the graph and unity. 
The value a = 1 is of course not a fractional order but, as we shall see, it is associated with a power- 
law tail of exponent minus one. 1 The minimum slope can be easily found. Indeed, F(f3) takes 
large values when (3 is large negative. In this case the main contribution to -F(/3)comes from the 
shortest interval of the partition with the length of the order of 7 2n . Hence, F(f3)/\@\ — > —2 In 7 
in the limit /3 — > —00. This gives the following lower bound of the differentiation order: 

Cn , 1-F" 1 (pln2) In 2 , . 

lim S = lim ( - '- = -— « 0.3777. (87) 

p— >oo p p^oo p 2 In 7 

We have already observed that, because the Cp graph does not pass through the origin, substantial 
values can be expected for the multifractality parameter C(a). The actual values of C, associated 
to values of p ranging from 1 to 3 by increments of 0.5 are shown in Tabled together with the 
number n of scaling octaves needed determined by nC ~ 10 (cf. Eq. (|86[)). In practice, on a 32 
bit machine, we are limited to about 25 octaves of dynamical range in resolution over the interval 
[— 1, 1]. This should be enough to observe power-law tails. 

Indeed, Fig. shows five instances of cumulative probabilities of fractional derivatives with 
power-law tails, corresponding to the values of the exponent p listed in Tabled The corresponding 
order of differentiation a ranges between 0.65 and 1. 2 Since the function u{x) which we are 
analyzing is not periodic, we resort to the Hann windowing technique employed previously in 
Ref. pQ (Section 13.4). Also, we use rank ordering to avoid binning. 

The power-law behavior observed is consistent with the phenomenological theory presented 
in Section 13.11 the residual discrepancies being due to the resolution of 2 25 bins. 

4 Concluding remarks 

We have found solid numerical evidence for the presence of power-law tails in the cumulative dis- 
tribution of fractional derivatives for the integral u{x) of the invariant measure of the Feigenbaum 
map. Furthermore the exponents measured are consistent with those predicted by phenomenolog- 
ical arguments from the spectrum of singularities. Since we have a fairly deep understanding of 

1 This is a consequence of £1 = 1 — F _1 (ln2) = 1. (Putting f3 = into the partition function Zk(f3) ~ 
E^o 1 lAfV, we have F k (0) = lnZ k (0)/* = ln2.) 

2 As predicted by the theory, when a is too small, e.g. for a — 0, no power-law tail is observed. 
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Figure 5: Cumulative probabilities of absolute values of fractional derivatives of various orders 
a = Cp*/P* f° r the Feigenbaum invariant measure. Each function displays a power-law tail 
with an exponent fairly close to the predicted value — p*. Insets: corresponding ( p graphs, (a) 
a = 1.0 (p* = 1.0). (b) a = 0.83 (p* = 1.5). (c) a = 0.74 (p* = 2.0). (d) a = 0.69 (p* = 2.5). (e) 
a = 0.65 (p* = 3.0). 
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the structure of the attractor, thanks in particular to the thermodynamic formalism, a reasonable 
goal may be to actually prove the results. The main difficulty is that the operation of fractional 
derivative is non-local. However, we believe that a rigorous analysis here is still possible due to 
the quite simple spectral structure of the dynamical system corresponding to the Feigenbaum 
attractor. 
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